Draft version February 5, 2008 

Preprint typeset using I^T^X style cmulatcapj v. 05/04/06 

@ 



The 10 /im infrared band of silicate dust: 
A laboratory study comparing the aerosol and KBr pellet techniques 

A. Tamanai and H. Mutschke 

Astrophysical Institute and University Observatory, Friedrich-Schiller-University Jena, Schillergasschen 3, D-07745 Jena, Germany 



o 
o 

(N 
Or 

00 



m 

(N 
ON 

o 
o 



J. Blum 

Institut fiir Geophysik und Extraterrcstrische Physik, Technische Universitat Braunschweig, Mendelssohnstrasse 3, D-38106 

Braunschweig, Germany 

AND 

G. Meeus 

Astrophysikalisches Institut Potsdam, An der Sternwarte 16, D-14482 Potsdam, Germany 
Draft version February 5, 2008 

ABSTRACT 

The profile of the silicate 10 /im IR band contains important information about the evolutional stage 
of dust in circumstcllar environments and the possible ongoing process of planetesimal formation. 
In order to extract this information, the observed band profiles are compared with calculated or 
laboratory-measured absorption cross sections of amorphous and crystalline grains with different sizes 
and compositions. We present in this study the first laboratory measurements of the 10 /urn band 
profiles of nonembedded, i.e. free-flying, particles of amorphous and crystalline Mg2Si04 (with two 
different particle shapes), amorphous and crystalline MgSiOa, and crystalline olivine. We compare the 
spectra with those measured on embedded grains and discuss the potential of the new experimental 
method for comparison with observed spectra, as well as for future studies of agglomeration and 
surface manipulation of the grains. 

Subject headings: circumstellar matter - stars: formation - stars: experiment - infrared - lines and 
bands - silicates 
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1. Introduction 

Dust particles are important players in astrophysical 
processes such as star and planet formation. Changes in 
size, composition, and crystallinity of the dust grains in- 
dicate dynamical processes (e.g. in accretion disks) and 
hold a key to understanding specifically the early for- 
mation of terrestrial (Earth-like) planets in the environ- 
ments of young stars. The most important technique, 
which makes it possible to obtain information about the 
dust populations in different objects (and nowadays even 
at different regions within one object), is infrared spec- 
troscopy. The most important role here is played by the 
10 /urn Si-0 stretching band of silicate dust because sili- 
cates are the major dust component in most dusty media 
in space, and this band can be observed from the ground. 

In fact, solid particles in the interstellar medium con- 
sist mainly of amorphous silicate dust grains (see the 
review bv IWhitteU Il992j) . whereas presence of crys- 
talline silicates has been confirmed via mid-infrared ob- 
servations in circumstellar disks around the relatively 
bright young Herbig A e/Be stars (e.g. IMalfait et al " 
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19981 IMeeus et al.l200lL for evolved stars (IMolster et al. 
200% . the T Tauri stars (TTSs: iMeeus et all 12003). 



the TTS Hen 3-600 (jHonda et al.l j20031. andjhe Vega- 
type stars HP 14526 3 and /3-Pic ie.g. lHonda et al]l2004: 
lOkamoto et al.ll2004|) . Analysis of the infrared emission 
bands and especially the 10 fj,m band profile of such sil- 
icate grains has allowed one to find indications of grain 
growth and crystallization that trace the early evolution 
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of young circumstellar disks toward planetary systems 
l)van Boekel et aLl l2""005). However, the analysis of obser- 
vational infrared spectra is very complex for an accurate 
interpretation, since not only mineralogical properties 
such as crystallinity and chemical composition determine 
the shape and intensity of the distinctive spectral feature, 
but also size, shape, and agglomeration of the grains ex- 
ert an influence. Thus, exhaustive laboratory studies of 
these influences on the silicate dust bands are essential 
for an appropriate interpretation of the observed spectra. 

Currently, band profiles that are used for comparison 
with observed spectra are mostly calculated, assuming 
simple geometrical models such as spherical or ellipsoidal 
grain shapes. The predictions of these calculations are 
uncertain because in reality the grain shape might be ir- 
regular. So far, laboratory measurements on real particle 
ensembles were not able to provide exact band profiles 
for a direct comparison due to the use of an embedding 
medium (KBr) in the sample preparation, which changes 
the band profile substantia lly by the influence of its elec- 
tromagnetic polarization (Fab ian et al.l 12001). In this 
Letter, we present a new method for laboratory mea- 
surements of extinction spectra avoiding the influence of 
electromagnetic interaction with solid embed ding media. 
Applying the aerosol technique (Hinds 1999), we obtain 
for the first time extinction spectra in the mid-infrared 
region of both amorphous and cr ystalline silicate grain s 
(for amorphous Si02 particles, see lTamanai et aD r2006). 
which are directly comparable to observed spectra. We 
investigate the differences to KBr-measured band profiles 
by deriving peak positions from the band profiles mea- 
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sured with both methods. Furthermore, we use electron- 
microscopic investigation of the grains to associate mor- 
phological (size, shape, agglomeration) influences on the 
measured profile. 

2. Experimental procedures and samples 

Regarding the aerosol measurements, we used a White- 
type long-path infrared cell (MARS-8L/20V, Gemini Sci- 
entific Instr.) of about 20 m path length, attached to 
a Fourier transformation infrared spectrometer (FTIR; 
Bruker ver. 113), to measure the extinction of the in- 
frared radiation by the silicate dust particles suspended 
in an aerosol cloud filling the cell (diameter about 15 cm, 
length 0.6 m). By multiple reflection between the two 
gold mirrors at both ends of the cell, the IR radiation 
passed 32 times through the particle cloud, thus en- 
hancing the sensitivity of the detection to optical depths 
around unity in the 10 /im band. An available dust flow 
generator (Palas RBG 1000), which utilized a rotating 
brush to disperse a pressed powder from a dust stor- 
age hole into a nitrogen gas stream, produced the dense 
aerosol that was led through stainless steel tubes toward 
the cell. In order to concentrate small-sized particles in 
the aerosol, we placed a two-stage impactor between the 
dust flow generator and the cell. The first stage retained 
clumps of particles above a certain size limit whereas 
the second stage densified the aerosol up to a concen- 
tration of 10 6 particles per cubic centimeter. A fraction 
of the aerosol was sampled through polyester membrane 
filters for scanning electron microscope (SEM) inspec- 
tion of the particle morphology. After filling the aerosol 
into the cell, we stopped the aerosol flow from the gen- 
erator and trapped the suspended particles in the cell 
during the spectroscopic analysis. Particles of 1 fim size 
and smaller are sufficiently coupled to the gas at atmo- 
spheric pressure to form a stable cloud on the timescale 
of hours. Larger particles and clumps of particles settled 
to the bottom of the cell. A small gas flow (50 liters 
hr_i) from both ends of the cell toward the outlet valve 
was maintained during the measurement in order to pro- 
tect the mirrors. This caused a continuous loss of aerosol 
and therefore a reduction of the particle column density 
during the experiment. Measurements shown here were 
usually taken after 10 minutes in order to allow for some 
sedimentation of remaining clumps. Apart from the de- 
creasing absorption strength, there were no changes of 
the band profiles observed after this time. 

In order to make a comparative study, we have ap- 
plied the classical KBr pellet technique where a granular 
sample (e.g. silicate particles) was mixed with potas- 
sium bromide powder (KBr), which has high transmis- 
sion throughout the mid-infrared range, in a mass ratio 
of 1:500. After very careful mixing in an agate mortar 
for segregating the particles as much as possible, 0.2 g of 
each mixture were pressed with 10 ton of force load to 
make about 0.5 mm thick clear pellets with a diameter 
of 1.3 mm, the transmission of which was measured in 
the FTIR spectrometer. 

We investigated six powdered silicate samples of five 
different sorts with respect to composition and crystal 
structure. The crystalline Mg2SiC>4 (forsterite) sam- 
ples were commercial products whereas for the amor- 
phous silicates we made use of original samples that 
were produced in our laboratory by the sol-gel technique 



llJager et al.l l2003|) and glass melting l|Dorschner et alJ 
11995ft . Table 1 gives a list of the samples we used, to- 
gether with their properties. We subjected the samples 
to a sedimentation process in a solvent (acetone), which 
allowed us to obtain a size fraction leq 1.0 /xm in diam- 
eter except CF2 and AF samples, which were originally 
small enough for the experiments. 

3. Results 
3.1. Aerosol vs. KBr spectra 

The extinction spectra for all samples are shown in Fig- 
ure 1 plotted as a normalized extinction efficiency versus 
the wavelength. The reason for normalizing all spectra to 
the maximum extinction values is that the aerosol mea- 
surements are not quantitative (see also § 4). 

From Figure 1, it is obvious that for the crystalline 
silicates, the aerosol spectroscopy reveals considerably 
different band positions and also bandwidths when com- 
pared to the KBr measurements. The band positions 
and their differences are given as inserted tables in the 
plots. Especially for the strong bands of the olivine-type 
crystals (CF1, CF2, CO) at -9.8 and ~11 fim, the spec- 
tra measured with the KBr pellet technique are shifted 
to longer wavelengths by amounts of up to 0.24 /xm. 
Smaller features which do not shift very much, change 
their appearance from isolated bands into shoulders and 
vice versa due to the shift of the strong bands. For crys- 
talline enstatite (CE), the shifts are not that strong (typ- 
ically 0.1 izm) but affect all the band components. For 
the amorphous materials (AF and AE), there is a clear 
detectable shift with AE, but a distinct shift is not seen 
with AF. 

The fact that strong bands are more severely influenced 
than weaker ones is expected. The strong bands are 
known to be due to geometrical resonances of the grains, 
which, for example, for spherical grains occur at e r 2e m 
where e and e m are the dielectric functions of the silicate 
and the embedding medium (1.0 for air, 2.3 for KBr), re- 
spectively. This can probably explain the detected shift 
toward longer wavelengths, where e decreases. For some 
theoretical simulations of band profil es for fo r sterite par- 
ticles in vacuum and KBr, see lFabian et alJ l)2001f >. On 
the other hand, the aggregation state of the particles in 
the aerosol and the KBr is probably different because 
of the different dispersion methods and the pressing of 
the pellets. Therefore, it may cause an additional influ- 
ence on the band profiles. Unfortunately, the agglomer- 
ation state in the KBr remains unknown, consequently 
the strength of this effect cannot be explored. 

Although in amorphous silicates negative values of e 
and therefore geometrical resonances do not occur, a 
small band shift is expected for particles in KBr com- 
pared to aerosol measurements. Calculations for the con- 
tinuous distribution of ellipsoids (CDE) model and with 
Mie theory give AA = 0.04 and 0.27 /im, respectively, for 
AE and AA = 0.02 and 0.20 /mi, respectively, for AF. It 
may have to do with the porous structure of the sol-gel 
material. Nevertheless, this shift is seen for AE but not 
for AF, which is currently not yet understood. 

3.2. Effect of grain shape 

Figure 2 shows transmission electron microscope 
(TEM) and SEM images of CF1 and CF2 samples. The 
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TABLE 1 



Properties of the silicate samples. 


Sample 


Chemical Formula 


Product Information 


Size (/^m) 


Shape 


CF1 


crystalline Mg2Si04 


Alfa Aesar Johnson 


d<1.0 


irregular 


CF2 


crystalline Mg2Si04 


Marusu Yuuyaku 


d^0.2 


ellipsoid 


AF 


amorphous Mg2.3Si04 


sol-gel process 


d<1.0 


irregular 


C'E 


crystalline MgSiOs 


melting 


d<1.0 


irregular 


AE 


amorphous MgSiOs 


melting/ quenching 


d<1.0 


irregular 


CO 


crystalline Mgi.9Feo.iSi04 


natural (San Carlos) 


d<1.0 


irregular 



9.86 


10.01 


0.15 


10.14 






10.43 


10.41 


0.02 


11.06 


11.23 


0.17 


11.88 


11.89 


0.01 




0.5 



Peaks: Cry. forsterite (irreg.: CF1) 
AS(nm) KBr(|.im) AX(|jm) 



Peaks: Cry. forsterite (ellip.:CF2) 
AS(jim) KBr(jim) AA(nm) 



9.81 


9.97 


0.16 


10.14 


10.41 


0.27 


10.83 


11.07 


0.24 


11.88 


11.88 







Peaks: Amo. MgzSiO< (AF) 
AS (nm) KBr (jim) M (nm) 
10,15 10.12 0.03 



0.5 



0.5 





Aerosol 




.... KBr 


/ / 




,7 Peaks: Amo. MgSiOs (AE) %N. 


// AS(nm) KBr(|.im) 
/ 9.6 9.78 
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A'k[\.im) 
0.18 


Aerosol *y t 

KBr / \; 


Peaks: Cry. enstatite (CE) 
AS(nm) KBr(nm) M(nm) 
9.23 9.21 0.02 
/>". 9.89 0.02 
.■/A\ 0.19 


,/ AS(jim) KBr(jim) 
/7 10.54 10.65 


A^(jim) ^ \ \ 
0.11 \\ 


-•/ 11.08 11.14 


0.06 V'.. 
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Fig. 1. — Comparison of the Si-O stretching band profiles of crystalline forsterite, amorphous Mg2Si04, crystalline enstatite, amorphous 
MgSiOs, and crystalline olivine obtained from the aerosol (AS) and KBr measurements. 



SEM pictures of the aerosol sampling demonstrate that 
most of the particles are in small and rather compact 
aggregates. The compactness indicates that they are 
clumps of the original powders and not agglomerates that 
have been grown in the aerosol. Although the clumps are 
less than 3 /im in size, the clumping should have an effect 
on the spectra. The close-packed clusters should tend to 
cause a broader band than the elongated ones as we have 
demonstrated by discrete dipole approximation (DDA) 
calculations for clusters of spherical SiC>2 particles (for 
more details see lTamanai et aLll200(Tfl . 

Although it is not possible to obtain the exact aggre- 
gate morphology, as the TEM images are not in direct 
view of the aerosol sampling, more detailed grain shapes 
can be examined. The images clearly show that CF2 



particles are rather roundish (ellipsoidal) , whereas those 
of CF1 are irregular with sharp edges. The band profiles 
measured for these samples clearly differ in the sense that 
the peaks of the geometrical resonances are shifted by up 
to 0.23 /im toward shorter wavelengths for the roundish 
CF2 grains. Although it cannot be directly proven that 
this is a pure consequence of the grain shape, it is plau- 
sible because the geometrical resonances of a sphere are 
situated at such shorter wavelengths (compare Fig. 3). 
In contrast, irregular grains generally seem to produce 
geometrical resonances at longer wavelengths, i.e. closer 
to the transverse optical lattice frequencies. 

3.3. Comparison with observations 
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(a) 



1.0|am 




Fig. 2.— (a,c) TEM and (b,d) SEM and images of the crystalline forsterite particles of samples (a,b) CF1 and (c,d) CF2. The SEM 
pictures show typical aggregates of 1-3 /im in size that have been directly sampled from the aerosol. The black dots in the SEM pictures 
are holes of the polyester filter having a diameter of 0.4 fim. The TEM pictures taken on the original powders give an impression of the 
grain shape. Note that sample CF2 has elliptical grains whereas those of sample CF1 are irregular. 

HD1 13766 differ in their morphological properties from 
the small roundish grains like CF2. 

4. Conclusion: comparison of the aerosol and KBr techniques 
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Fig. 3. — Top: Comparison of the HD113766 emission feature 
with the CF1 band profile from the aerosol and KBr measurements. 
Bottom: Calculated band profiles for forsterite using the CDE and 
DHS models. The light dotted lines denote the peak positions of 
the aerosol result. 



Finally, we demonstrate the importance of our new 
data by applying them to an observed spectrum. Figure 
3 shows the emission spectrum o f the crystalline sil icate- 
rich Vega-type star HD113766 <|Schutz et alJl2005|) . to- 
gether with the band profiles of CF1 measured via 
both the aerosol and the KBr measurements, as well 
as calculated band profiles of crystalline forsterite ( com- 
plex refractive in dex bv iServoin fc Pirioul Il973ft for 
the CDE model l|Bohren & Hu ffman 1983), and for 
the distribution of ho llow spheres (DHS) model (see 
Ivan Boekel et~alll2005j) . 

Evidently, the calculated band profiles cannot re- 
produce both the observed and the laboratory-measured 
data. The CDE model leads to band profiles that are 
far broader than the laboratory measurements. The 
DHS model predicts the structure of the observed band 
correctly but produces peaks at excessively long wave- 
lengths, particularly at 10.1 fim (similar to the CDE). 
In contrast to this, the aerosol measurements of the CF1 
sample reproduce quite well the positions of the main 
peaks in the observed spectrum at 9.9 and 11.0 fim. The 
peak at 11.3 fim, however, is not reproduced. This may 
indicate larger grains in the disk of the Vega-type star. 
We should note here that the measured spectrum of the 
natural olivine grains (sample CO) provides a similarly 
nice match whereas the band profile of sample CF2 fails 
to reproduce the observed features (see Fig. 1). This 
may indicate that the forsterite- or olivine-type grains in 



We have applied the aerosol technique to IR spectro- 
scopic measurements of the 10 fim band profile of crys- 
talline and amorphous silicates and demonstrated its ap- 
plication to obtaining realistic band profiles that can be 
compared directly to observed spectra for the first time. 

Compared to the traditional KBr technique, the ex- 
tinction measurement in aerosol has two disadvantages: 
it requires about 100 mg of material and it is, unfor- 
tunately, not quantitative. The column density of the 
dust grains along the spectrometer beam can hardly be 
estimated or measured to an accuracy higher than a fac- 
tor of 2. A better accuracy could be achieved only by 
accompanying quantitative KBr measurements or by a 
determination of the material's optical constants, for ex- 
ample, by reflection measurements. In such a case and 
given that a theoretical model can be adjusted to repre- 
sent correctly the morphological properties of the partic- 
ulate, the aerosol spectra can be renormalized. 

Conversely, there are distinct advantages of the aerosol 
method. First, it lacks the problem of the embedding 
medium, which is inherent to the KBr method. A sec- 
ond major advantage is that the structures of the grains 
can be analyzed by either in situ or ex situ microscopic 
imaging. This allows investigation in detail of the in- 
fluence of the grain morphology on the spectra and will 
hopefully help to calibrate theoretical approaches. 

Concerning practical aspects, the preparation of the 
aerosol is easier than that of a KBr pellet with an ideal 
homogeneous grain distribution. Large clumps of parti- 
cles can be receded or settled down within a short time. 
Even separation of different grain sizes is uncomplicated 
with the impactor technique. An important point is that 
collision timescales in the aerosol can probably be made 
smaller than the duration of an experiment, so that frac- 
tal aggregate formation could be directly studied by IR 
spectroscopy. So far, the particle densities were too small 
in our experiment to observe grain growth by aggrega- 
tion. With the KBr technique, this is hardly possible, 
because agglomerates are modified or destroyed by the 
preparation process. Finally, with the aerosol technique, 
it should be possible to manipulate the dust in situ by 
surface condensation or irradiation, so that its influences 
on the IR spectra can be measured well. 
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